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SUPERCRITICAL PRESSURE REGULATION OF VAPOR COMPRESSION 
SYSTEM BY REGULATION OF ADAPTIVE CONTROL 

BACKGROUND OF THE INVENTION 
[1] The present invention relates generally to a system for regulating the high 

pressure component of a transcritical vapor compression system by employing adaptive 
control. 

[2] Traditional vapor compression systems are designed to run at zero superheat, or 
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subcritically are commonly optimized by choosing load characteristics (high side heat 
rejection), input characteristics (low side heat absorption), and refrigerant characteristics 
(high/low side refrigerant temperature and superheat). 

[3] Chlorine containing refrigerants have been phased out in most of the world as 

they have a negative effect on ozone. Hydrofluoro carbons (HFCs) have been used as 
replacement refrigerants, but these refrigerants may also have negative effects. "Natural" 
refrigerants, such as carbon dioxide and propane, have been proposed as replacement 
fluids. Carbon dioxide has a low critical point, which causes most air conditioning 
systems utilizing carbon dioxide to run partially above the critical point, or to run 
transcritical, under most conditions. The pressure of any subcritical fluid is a function of 
temperature under saturated conditions (when both liquid and vapor are present). 
However, when the temperature of the fluid is higher than the critical temperature 
(supercritical), the pressure becomes a function of the density of the fluid. 

[4] When a vapor compression system is run transcritical, it is advantageous to 

regulate the high pressure component of the system. By regulating the high pressure of 
the system, the capacity and/or efficiency of the system can be controlled and optimized. 
By operating the system transcritically, the pressure and the temperature of heat rejection 
can both be independently controlled. 

[5] Adaptive control can be employed to modify variable coefficients in an adaptive 

control algorithm. By modifying the variable coefficients in the adaptive control 
algorithm, the optimal variable setpoint reference that achieves maximum capacity can be 
obtained. 



Express Mail No. EV221417044US 



60246-189 
10580 



In a prior vapor compression system, receding horizon control using Recursive 
Least Squares (RLS) has been used for model identification. Systems have also been 
optimized by adjusting control variables and solving the system through direct matrix 
inversion. Both of these approaches automatically adjust the variable speed of the heat 
pump, the blower speed, and the evaporator superheat of a vapor compression system to 
optimize the coefficient of performance of the system. 

SUMMARY OF THE INVENTION 

A vapor compression system in basic form includes a compressor, a gas cooler, an 
expansion device, and an evaporator. Refrigerant is circulated though the closed circuit 
cycle. Carbon dioxide is used as the refrigerant. As carbon dioxide has a low critical 
point, systems utilizing carbon dioxide as a refrigerant usually require the vapor 
compression system to run transcritical. When the system is run transcritical, it is 
advantageous to regulate the high pressure component of the vapor compression system 
to control and optimize the capacity and/or efficiency of the system. The overall 
efficiency of the system is determined by comparing the amount of useful energy 
extracted by the heat rejecting heat exchanger to the amount of energy expended to 
compress the refrigerant and to run any ancillary components of the system, such as heat 
exchanger fans or pumps. 

As the environment and the system changes over time, the high pressure that 
provides the maximum coefficient of performance of the system changes. Adaptive 
control is employed to modify the model that operates the system to continually optimize 
the coefficient of performance. The model is determined by an adaptive control 
algorithm including variable coefficients. As the system changes over time, the model 
that operates the system is modified to optimize the coefficient of performance. As the 
variables of the adaptive control algorithm change, the model changes. A control of the 
heat rejecting heat exchanger is then adjusted based on the modifications to regulate the 
high pressure of the system and therefore the coefficient of performance. 

In a first example, Least Mean Squares (LMS) is used to modify the variables of 
the adaptive control algorithm to optimize the coefficient of performance. In the first 
step, a system identification error is computed using a gradient descent methodology. In 
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the second step, the model is adapted using the system identification error information. 
The model is modified such that the output of the model is substantially equal to the 
output of the system, reducing the system identification error to zero. A control is 
adjusted based on the adaptive control update to adjust the high pressure of the system in 
the heat rejecting heat exchanger to obtain the maximum coefficient of performance. 
[10] In a second example, the coefficient of performance is optimized by employing a 

slowly varying periodic excitation method to seek extreme conditions. Intelligent 
excitation and signal manipulation and filtering are employed to achieve an adaptive 
control update. 

[11] A third example employs triangularization to find the optimal coefficient of 

performance of the system. A triangle having three setpoints is established: one point is 
on the left hand side and has a positive slope, one point in on the right hand side and has 
a negative slope, and a midpoint lies between these points. After first establishing a 
triangle including three points, the adaptive control algorithm constricts the control 
variables and focuses on the old midpoint to define a new middle point. 
Triangularization is repeated until the maximum coefficient of performance is obtained. 
The system is then run at this input to maximize the coefficient of performance. 

[12] These and other features of the present invention will be best understood from the 

following specification and drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[13] The various features and advantages of the invention will become apparent to 

those skilled in the art from the following detailed description of the currently preferred 
embodiment. The drawings that accompany the detailed description can be briefly 
described as follows: 

[14] Figure 1 illustrates a schematic diagram of a prior art vapor compression system; 

[15] Figure 2 illustrates a thermodynamic diagram of a transcritical vapor compression 

system; 

[16] Figure 3 schematically illustrates an adaptive control system employed to 

maximize the coefficient of performance of a transcritical vapor compression system 
employing the Least Mean Squares Approach; 
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[17] Figure 4 schematically illustrates a system employed to maximize the coefficient 

of performance of a transcritical vapor compression system employing slowly varying 
periodic excitation; 

[18] Figure 5 schematically illustrates a system employed to maximize the coefficient 

of performance of a vapor compression system employing triangularization; 
[19] Figure 6 schematically illustrates the initialization of the system of Figure 5; and 

[20] Figure 7 schematically illustrates the triangularization found of the system of 

Figure 5. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[21] Figure 1 illustrates a vapor compression system 20 including a compressor 22, a 

gas cooler/condenser 24, an expansion device 26, and an evaporator. Refrigerant 
circulates though the closed circuit system 20. 

[22] The refrigerant exits the compressor 22 at high pressure and high enthalpy and 

flows through the gas cooler 24. In the gas cooler 24, the refrigerant loses heat, exiting 
the gas cooler 24 at low enthalpy and high pressure. A fluid medium, such as water or 
air, flows through a heat sink 32 and exchanges heat with the refrigerant passing through 
the gas cooler 24. For example, cooled water enters the heat sink 32 at the heat sink inlet 
34 and flows in a direction opposite to the direction of flow of the refrigerant. After 
exchanging heat with the refrigerant, the heated water exits through the heat sink outlet 
36. The refrigerant then passes through the expansion device 26, and the pressure drops. 
After expansion, the refrigerant flows through the evaporator 28 and exits at a high 
enthalpy and low pressure. 

[23] Preferably, carbon dioxide is used as the refrigerant. While carbon dioxide is 

illustrated, other refrigerants may be used. Because carbon dioxide has a low critical 
point, systems utilizing carbon dioxide as a refrigerant usually require the vapor 
compression system 20 to run transcritical. 

[24] The refrigerant exits the compressor 22 at high pressure and high enthalpy, shown 

by point A in Figure 2. As the refrigerant flows through the gas cooler 24 at high 
pressure, it loses heat and enthalpy, exiting the gas cooler 24 with low enthalpy and high 
pressure, indicated as point B. As the refrigerant passes through the expansion device 26, 
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the pressure drops, shown by point C. After expansion, the refrigerant passes through the 
evaporator 28 and exits at a high enthalpy and low pressure, represented by point D. 
After the refrigerant passes through the compressor 22, it is again at high pressure and 
high enthalpy, completing the cycle, 
[25] The overall efficiency of the system 20 is determined by comparing the amount of 

useful energy extracted by the gas cooler 24 to the amount of energy expended to 
compress the refrigerant in the compressor 22 and run any ancillary components of the 
system 20, such as heat exchanger fans or pumps. Therefore, the system performance 
(z(u)), or coefficient of performance, is defined as: 

z(u) = L Useful Energv 

S Parasitic Energy (Equation 1) 

[26] The variable u is the pressure setpoint which is controlled by the control input to 

the expansion device 26. Therefore, the coefficient of performance of the system 20 is a 
function of the pressure setpoint. 

[27] When a vapor compression system 20 is run transcritical, it is advantageous to 

regulate the high pressure side of the system 20. By regulating the high pressure of the 
system 20, the coefficient of performance of the system 20 can be controlled and 
optimized. By operating the system 20 transcritically, the pressure and the temperature 
of heat rejection can be independently controlled. 

[28] Adaptive control of a transcritical vapor compression system 20 offers variable 

setpoint references and variable algorithm parameters that are used to achieve the 
setpoints. The variable setpoint references modify themselves to facilitate optimal 
operation of the system 20 regardless of external influences. That is, the setpoint 
references are modified to obtain the optimal coefficient of performance. By varying the 
parameters of the control algorithm, the impact of external influences on the coefficient 
of performance to achieve the setpoints is minimized. 

[29] The system 20 is operated by a model. As the environment and the system 20 

change over time, the high pressure which provides the maximum coefficient of 
performance of the system 20 changes. If the model is not updated as the system 20 
changes, the model does not property operate the changed system 20, resulting in a non- 
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Optimal coefficient of performance. In the present invention, adaptive control is 
employed to modify the model to continually optimize the coefficient of performance of 
the system 20. The model is determined by an adaptive control algorithm including 
variable coefficients. As the system 20 changes over time, the variable coefficients of the 
adaptive control algorithm are changed to modify the model to optimize the coefficient of 
performance. The variables are modified based on influences such as ambient conditions 
and user-selected loading on the system 20. 

[30] Most adaptive control approaches include two steps: 1) initialization at a setpoint 

near to the expected optimum and 2) subsequent excitation to determine which direction 
in parameter space further maximizes the set of variables in question. 

[31] Preferably, there are two independent control loops in the system 20. A PID loop 

for setpoint control in the gas cooler 24 maintains regulatory control of the high side 
outlet temperature. An LMS loop for optimization of the coefficient of performance 
changes the setpoint according to how the system 20 is running. Both loops are coupled 
by the heat pump. Coupling is minimized by operating the loops on different time scales. 
The PID loop runs relatively quickly and the LMS loop runs relatively slowly. 
Therefore, the LMS loop does not respond to changes in the high side temperature input 
or output. The excitation of the system 20 for system identification runs between the 
speed of the PDD loop and the speed of the LMS loop. 

[32] In a first example, as shown in Figure 3, Least Mean Squares (LMS) is used to 

modify the variables of the adaptive control algorithm to optimize the coefficient of 
performance. In the first step, a system identification error is computed using a gradient 
descent methodology. The system identification error is the difference between the 
output of the system and the output of the model. In the second step, the model is 
adapted using the system identification error information. 

[33] The system is identified by a Least Mean Squares (LMS) approach. The system 

20 is first excited (du) by adjusting the control 38 (shown in Figure 1). The output of the 
system (dy) is then compared to the output of the system model (dy'). The output of the 
system (dy) is sensed by a sensor 40. The system model (O) is the slope of the 
coefficient of performance (z(u)) with respect to the pressure setpoint (u). Therefore, the 
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output of the system model (dy*) is the product of the system model (O) and the first 
excitation (du). 



which further maximizes the coefficient of performance (z(u)). The excitation (du) must 
be continuous to account for slow changes in the coefficient of performance. The 
frequency of the excitation (du) depends on system dynamics. A static relation is chosen 
for the model identification, and the response of the system 20 must appear to be static 
with respect to the control signal. In one example, the excitation (du) is a low amplitude 
sinusoidal input. The amplitude is chosen so that the excitation (du) perturbs the 
coefficient of performance (z(u)) outside levels of instrumentation noise. 

[35] The model (O) is determined by an adaptive control algorithm having variable 

coefficients. The system 20 and the optimal coefficient of performance (z(u)) of the 
system 20 change over time. The variable coefficients of the adaptive control algorithm 
are modified based on the state of the system 20 to adapt the system 20 to optimize the 
coefficient of performance. The variable coefficients are modified based on influences 
such as ambient conditions and user-selected loading on the system. The Least Mean 
Squares approach is employed to adapt the system model (O) so the output of the model 
(dy') matches the output of the system (dy). 

[36] After excitation (du) of the system, the output of the system (dy) is then compared 

to the output of the system model (dy') to determine the system identification error (e). 
The system identification error (e) is the difference between the output of the system (dy) 
and the output of the model (dy'), or 



[34] 



The system 20 is excited bi-directionally to determine the direction in space 



e = dy-dy' = dy-Odu. 



(Equation 2) 



[37] 



The system identification error (e) is then squared to obtain a positive quadratic 
function of model accuracy with respect to parameters of the model: 



E[e^] = E[dy'^] + <D^[dudu'^]0 - 20E[dydu]. 



(Equation 3) 
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[38] By differentiating, the gradient of the quadratic relation is obtained: 

V= gEre^l = 20E[dudu'^] - 2E[dydu]. (Equation 4) 

SO 

[39] By forcing discrete-time model updates in the direction of the gradient, the system 

identification error (e) is minimized (steepest descent methodology), 

[40] A model update parameter (fi) is employed to generate the identified model 

update. The model update parameter (]I) affects stability and convergence time of the 
system identification. The model update parameter (p) is dependent on the excitation of 
the system (du) and is defined as: 

p; = \i (Equation 5) 

e+du^ 

[41] The user defined adaptation constant (|a) is typically chosen small, for example 

between 0.0 and 1.0. The normalization parameter (s) alleviates numerical difficulties 
when the model update parameter (]I) is small and chosen based on expected minimum 
levels of the control signal. Normalization reduces noises and unknown system dynamics 
that occur during high amplitude excitation which can affect the model update parameter 

(?)• 

[42] The identified model update (<Dk+i) is defined as: 

Ok+1 = Ok + F (-Vk) (Equation 6) 

[43] where k is a point in time. That is, the update model (Ok+1) is the sum of the 

current system model (Or) and the product of the model update parameter (p) and the 

negative gradient (-Vk). 
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[44] After identification of the response to an input, the optimal adaptive control 

update to the input (du) to the system 20 is derived using the optimal control theory. The 
known quadratic cost function is defined as: 

J(u) = iJqz + du'^adu. (Equation 7) 

[45] In the quadratic cost function, q and a are user defined compensation weights and 

z is the implied system performance parameter of the vapor compression system 20 near 
the mean operating point. 

The system performance parameter (z) is defined as: 

z = zo + Odu. (Equation 8) 

[46] In the system performance parameter equation, zo is the mean operating point of 

the system performance parameter (z). The quadratic cost function J(u) is employed to 
penalize both the system performance parameter (z) and a control derivative. The 
parameters of the cost function J(u) are chosen while balancing the rate of convergence to 
the optimum and stability. 

[47] By solving for the extremum of the quadratic cost function J(u), the optimal 

adaptive control update is obtained: 

du = - Oq (Equation 9) 

O^qO + a. 

[48] By employing the identified model update Ok+i in the optimal adaptive control 

update equation, the optimal adaptive control update can be obtained. The optimal 
adaptive control update (du) adapts the model (O) such that the output of the model (dy') 
is substantially equal to the output of the system (dy), reducing the system identification 
error (e) to zero. The control 38 is adjusted based on the adaptive control update (du) to 
adjust the high pressure of the system 20 in the gas cooler 24 to obtain the maximum 
coefficient of performance. In one example, the control 38 adjusts a variable speed pump 
for water flow through the gas cooler 24. 

[49] In a second example, shown in Figure 4, the coefficient of performance of the 

system 20 is optimized by employing a slowly varying periodic excitation method. 
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Intelligent excitation and signal manipulation, including signal filtering, are employed to 
achieve an adaptive control update in the direction of the extreme. 

[50] The approach begins with a sinusoidal excitation of the system 20 superimposed 

on the mean setpoint. The frequency of the excitation is slow compared to the system 
dynamics. The sinusoidal excitation is defined as: 

u' = u + Asin(cot) (Equation 10) 

[51] where u' is the signal input to the system 20 and u is the mean signal input to the 

system 20. 

[52] The system response to the sinusoidal excitation (u') is: 

z' = zo + z sin(cot) (Equation 11) 

[53] The response (z') is then filtered by a first dynamic filter Hi(s) (s is the Laplace 

variable). The first dynamic filter Hi(s) is high-pass in nature and removes the DC offset 
from the response (z'), resulting in a harmonic response. Therefore, the response (z') 
must have zero phase at the frequency of the sinusoidal excitation signal (u'). If the 
sinusoidal excitation signal (u') is chosen within the bandwidth of the system 20, the first 
dynamic filter Hi(s) can partially account for dynamic interaction of the system 20. The 
first dynamic filter Hi(s) also reduces instrumentation noise. 

[54] Assuming slow excitation, the harmonic response contains either an in-phase or 

out-of-phase sinusoidal response with respect to the sinusoidal excitation signal (u'). 
This phase characteristic facilitates the determination of the setpoint update in the 
direction of the extreme. If the system 20 is in phase, u is less than Uopt, where Uopt is the 
optimum solution; alternately, if the system 20 is out of phase, then u is greater than Uopt- 
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The harmonic response is then multiplied by the initial excitation u to demodulate 
the harmonic response. By demodulating the harmonic response with the initial 
excitation u, a relative static gain and the high-frequency harmonic information is 
obtained: 

Z*u = V2va: - !/2Uz'sin(2(0t) (Equation 12) 

The high frequency content is then filtered using a second dynamic filter HaCs), 
The second dynamic filter H2(s) is low-pass in nature and removes the oscillating factor 
\}/2nz' sin(2(ot)] to leave only the static gain (l/2uz'). The static gain is similar to the slope 
of the system 20 coefficient of performance curve z(u). An integral-type controller 
coupled to the second dynamic filter H2(s) drives the relative static gain (or slope) of the 
system response (z') to zero. The second dynamic filter H2(s) supplies a new mean 
control signal sufficient to drive the slope to zero and facilitates the stability and 
performance of the control algorithm. The steady state gain of the second dynamic filter 
H2(s) must match the amplitude of the excitation for stability. In addition, the pole of the 
second dynamic filter H2(s) is chosen to distinguish needed information from the 
demodulation process. The frequency response of the second dynamic filter H2(s) can be 
shaped to decrease the bandwidth of the closed-loop system, increasing the convergence 
over time. 

The filtered signal from the second dynamic filter H2(s) is then used as the input 
of the system 20. As the coefficient of performance z(u) is a function of the input (u') to 
the system 20, the coefficient of performance is modified to an optimal value. The 
control 38 is adjusted to adjust the high pressure of the system 20 in the gas cooler 24 to 
obtain the maximum coefficient of performance. The advantage of this method is that the 
system 20 identification does not interrupt the control, resulting in a better rate of 
convergance. 

Referring to Figure 5, a third example of the present invention employs 
triangularization to determine the maximum coefficient of performance (z(u)) of the 
system 20. Triangularization includes 1) discrete excitation of the system 20 and 2) 
memory of the system response to adaptively map the system 20. By intelligent 
movement of the control setpoint, a triangle can be established. 
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[59] Figure 5 illustrates a plot of the coefficient of performance z(u) as a function of 

the input (u) of the mean control variable of the system 20. A triangle having three 
setpoint conditions is established. One point z(ulhs) is on the left hand side of the graph 
and is a positive value because the curve is a positive slope, one point z(urhs) is on the 
right hand side of the graph and is a negative value because the curve has a negative 
slope, and a point z(umax) which indicates the maximum coefficient of performance z(u) 
lies between these points. Therefore, 

dz\ >0; and5z| <0; ->Umax e (ulhs, urhs) (Equation 13) 

u=uRHS 

[60] The history of the response of the system 20 is maintained to always retain a 

triangle of conditions having the extreme lying within. After first establishing a triangle 
including points z(ulhs), z(urhs), z(umid), [umid is a point between z(Ulhs) and z(Urhs)], 
the algorithm constricts the control variables and focuses on Umid to define a new middle 
point. For example, as shown in Figure 6, the system 20 is first brought to operating 
point ui. Movement is made from the first setpoint uj to a second setpoint U2. If this 
movement results in a decreased coefficient of performance, movement is made in an 
opposite direction. For example, as shown in Figure 6, movement to the left offers a 
decreased coefficient of performance. Movement is then made to the right which offers 
an increased coefficient of performance. A third set U3 point is chosen in a direction 
similar to the direction of movement to the second setpoint U2. This movement offers an 
increased coefficient of performance and the third setpoint U3 is chosen. However, if this 
movement offered a decreased coefficient of performance, the third setpoint U3 is chosen 
in an opposing direction. The identification of the triangle is accelerated if the proper 
choice of setsize is established. Therefore, if 

z(u2) > z(ui) ^ correct direction U3 > U2;; and if (Equation 14) 

z(u2) < z(ui) — ► wrong direction U3 < ui. 

[61] The initial stepsize must be chosen small so an initial step in the wrong direction 

will not degrade system 20 performance. The minimum stepsize is the minimum amount 
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of control variation needed to vary the coefficient of performance. The time in which 
each step is made is based on characteristics of the dynamic response of the system 20. 
However, a minimal stepsize must be designated to ensure that if the nonlinear 
coefficient of performance z(u) changes over time, the changes will be tracked by the 
algorithm (persistent excitation). The stepsize must also be limited to a maximum for the 
initialization process to limit it from overshoot from one side of the maximum to the 
other. 

[62] The choice of the stepsize is established by multiplying a gain, such as 

Fibonacchi's golden ratio (k=1.618), by the last successful stepsize: 



Aunew = kAuoid, k>l. (Equation 15) 

[63] As shown in Figure 7, initialization is continued until a triangle having a 

maximum performance u^-i situated within the setpoints Uk and Uk.2 is obtained. To verify 
this, the perturbations are monitored until a set of three setpoint conditions having one 
maximal value situated between two smaller values is obtained. When this condition is 
met, the following is true: 



z(Uk-i) > z(Uk.2); and z(Uk) < z(uk-i) triangle established. (Equation 16) 

[64] Once a triangle is established, initialization is complete. The algorithm then 

constricts the size of the triangle by focusing on the current midpoint Uk-i. New setpoints 
are chosen as logical perturbations from the midpoint Uk-i, while considering the farthest 
previous point from the midpoint Uk-i of the previously accepted triangle. By considering 
the farthest previous points, the triangle is reduced in size. Additionally, if the stepsize is 
chosen based on the previous stepsize, the algorithm deccelerates when nearing the 
maximum, facilitating stability. 

[65] The triangular restriction is repeated until the input (u) that obtains the maximum 

coefficient of performance z(u) is obtained. When the input (u) is obtained, the system 
20 is then run at this input (u) to maximize the coefficient of performance z(u) of the 
system 20. The control 38 is adjusted to adjust the high pressure of the system 20 in the 
gas cooler 24 to obtain the maximum coefficient of performance. 
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[66] Alternately, a single multiple input, multiple output (MIMO) adaptive controller 

can be utilized rather than the two independent controls PID and LMS loops to provide 
faster coefficient of performance convergence and adaptation of the water temperature 
regulation. Simultaneously regulation of the high side outlet temperature and 
maximization of the coefficient of performance is accomplished using a multi-variable 
form of the gradient descent control and adaptation utilized to optimize the coefficient of 
performance discussed above. The water outlet temperature is not invariant with respect 
to the expansion valve 26, and the coefficient of performance is influenced by the rate of 
energy extraction in the heat rejecting heat exchanger 24. 

[67] The foregoing description is only exemplary of the principles of the invention. 

Many modifications and variations of the present invention are possible in light of the 
above teachings. The preferred embodiments of this invention have been disclosed, 
however, so that one of ordinary skill in the art would recognize that certain 
modifications would come within the scope of this invention. It is, therefore, to be 
understood that within the scope of the appended claims, the invention may be practiced 
otherwise than as specifically described. For that reason the following claims should be 
studied to determine the true scope and content of this invention. 
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